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The fully self-consistent Random Phase Approximation (RPA) is constructed within 
the Unitary Correlation Operator Method (UCOM), which describes the dominant 
interaction-induced short-range central and tensor correlations by a unitary transfor- 
mation. Based on the correlated Argonne V18 interaction, the RPA is employed in 
studies of multipole response in closed-shell nuclei across the nuclide chart. The UCOM- 
RPA results in a collective character of giant resonances, and it describes rather well 
the properties of isoscalar giant monopole resonances. However, the excitation ener- 
gies of isovector giant dipole resonances and isoscalar giant quadrupole resonances are 
overestimated due to the missing long-range correlations and three-body contributions. 



1. Introduction 

One of the unresolved problems in the theory of nuclear structure is the descrip- 
tion of ground state properties and excitation phenomena for heavier nuclei, based 
on realistic nucleon-nucleon (NN) interactions which reproduce the NN scattering 
data 1 ' 2 - 3 . The use of these realistic interactions for solving the nuclear many-body 
problem is a challenging task. Presently, only light nuclei can be treated within ab 
initio schemes like Green's function Monte Carlo 4 , no-core shell model 5 , and cou- 
pled cluster method 6 . The Unitary Correlation Operator Method (UCOM), which 
describes the dominant short-range and tensor correlations explicitly by means of 
a unitary transformation 7 ' 8 > 9 > 10 ; allows for the use of realistic NN interactions in 
traditional nuclear structure methods. In contrast to other methods using unitary 
transformations, e.g. the unitary model operator approach n > 12 ; the correlation 
operators are given explicitly allowing for the derivation of a system-independent 
effective interaction operator VjjCOM- Although different by its construction, the 
correlated NN interaction Vucom is similar to the Vi ow -k low-momentum interac- 
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tion 13 . Both approaches lead to the separation of momentum scales, providing a 
phase-shift equivalent NN interaction in the low-momentum regime. 

Very recently, the correlated realistic NN interaction constructed within the 
UCOM framework, has been employed in Hartree-Fock (HF) calculations across 
the nuclide chart 14 . Based on the UCOM-HF ground state, we have constructed 
the random-phase approximation (RPA) for the description of low-amplitude col- 
lective excitations in atomic nuclei using correlated realistic NN interactions 15 . 
Various phenomenological RPA and quasiparticle RPA models have been very suc- 
cessful in the past, not only in studies of giant resonances and low- lying states 
(e.g. Rcfs. i 7 * 18 * 19 * 20 . 21 ^ but a i so m description of exotic nuclear structure of col- 
lective excitations in nuclei away from the valley of /3-stability 20 > 22 > 23 > 24 > 25 > 26 < 27 > 28 . 
The present study, however, provides the first insight into the collective excitation 
phenomena in closed-shell nuclear systems, based on the correlated realistic NN 
interactions. 



2. Unitary Correlation Operator Method (UCOM) 

The essential ingredient of the UCOM approach is the explicit treatment of 
the interaction-induced correlations, i.e. short-range central and tensor correla- 
tions 7 ' 8 ' 9 . The relevant correlations are imprinted into an uncorrelated many-body 
state I^F) through a state-independent unitary transformation defined by the uni- 
tary operator C, resulting in a correlated state |^) = C \ *f?) . An equivalent, tech- 
nically more advantageous approach, is based on using the correlated operators 
O = C^OC with uncorrelated many-body states. The short-range central corre- 
lations are described by a distance-dependent shift, pushing two nucleons apart 
from each other if they are closer than the core distance. The application of the 
correlation operator in two-body space corresponds to a norm conserving coordi- 
nate transformation with respect to the relative coordinate. This transformation is 
parameterized in terms of correlation functions for each (S, T) channel which are 
determined by an energy minimization in the two-body system. For purely repul- 
sive channels an additional constraint on the range of the central correlator is used 
(Ir _0 ' T_0 ' ) =0.1 fm 4 , cf. Rcf. 10 ). The details of the determination and parame- 
terization of the standard correlator are given in Ref. 10 . The tensor correlations 
between two nucleons are generated by a tangential shift depending on the spin 
orientation 8 . The size and the radial dependence is given by a tensor correlation 
function $(r) for each of the two 5 = 1 channels, whose parameters are also deter- 
mined from an energy minimization in the two-body system 10 . The range of the 
tensor correlation function is restricted through a constraint on the range measure 
1$ = J drr 2 d{r). If one would use for the description of finite nuclei the long- 
range tensor correlator that is optimal for the deuteron, an effective screening due 
to other nucleons would emerge through higher-order contributions of the cluster 
expansion 8 . In practical calculations based on two-body approximation, this prob- 
lem is resolved by restricting the range of the tensor correlation function, which 
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provides an effective inclusion of the screening effect without explicit evaluation of 
higher terms in the cluster expansion 9 . Recent studies within the exact no-core 
shell model 10 , show that I# ~ ' = 0.09 fm 3 leads to an optimal tensor corre- 
lator for the description of binding energies of 3 H and 4 Hc. In the present work we 
vary the range of the tensor correlator, i^ S_1,T_0 ^ = 0.07, 0.08, and 0.09 fm 3 , in 
order to probe its impact on the description of the global properties of collective 
excitation phenomena in atomic nuclei. The contributions of the tensor correla- 
tor in (S,T) = (1, 1) channel are one order of magnitude smaller 9 , and therefore 
neglected in the present study. 



3. Random-phase approximation in the UCOM framework 

Starting from the uncorrelated Hamiltonian for the A-body system consisting of the 
kinetic energy operator T and a two-body potential V, the formalism of the unitary 
correlation operator method is employed to generate the correlated Hamiltonian. By 
combining the central and tensor correlation operators, the correlated many-body 
Hamiltonian in two-body approximation is given by, 

Hucom = TW + + yPl = T + Vucom, (1) 

where the one-body contribution comes only from the uncorrelated kinetic energy 
TW = T. Two-body contributions arise from the correlated kinetic energy T^l and 
the correlated potential V^ 2 \ which together constitute the correlated interaction 
Kjcom 10 ■ More details about the evaluation of the two-body matrix elements 
for Vucom are available in Ref. 10 . Assuming spherical symmetry, the correlated 
realistic NN interaction is employed to solve the HF equations, i.e. to evaluate 
the single-particle wave functions and energies 14 . The UCOM-HF single-particle 
spectra are used as a basis for the construction of the p — h configuration space 
for the RPA. The RPA equations are derived from the equation of motion method 
using the quasiboson approximation 17 , 



A J B J \ fX v ' JM \ /l \ fX^ JM 



B *J A *J I 1 yv,JM 0-1 Y»' JM 



(2) 



where the eigenvalues uj v correspond to RPA excitation energies. The residual 
particle-hole interaction in A and B matrices includes the correlated realistic NN 
interaction Vucom in a fully consistent way with the Hartree-Fock equations. In 
addition, the multipole transition operators are consistently transformed by em- 
ploying the same unitary transformation as for the Hamiltonian. 

However, it turns out that the effect of the UCOM transformation of the transi- 
tion operators for monopole and quadrupole modes is negligible 15 . It is interesting 
to note, that the UCOM-RPA results are in agreement with the study of effective 
operators in the no-core shell model within the 2MI model space, where the B(E2) 
values are very similar for the bare and the effective operator which includes the 
two-body contributions 16 . An essential property of the present UCOM-RPA scheme 
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is that it is fully self-consistent, i.e. the same correlated realistic NN interaction is 
used in the HF equations that determine the single-particle basis, as well as the 
RPA residual interaction, and the multipolc transition operators are transformed 
consistently with Vucom • This means that the same unitary transformation of the 
realistic NN interaction, i.e. central and tensor correlation functions with the same 
set of parameters are systematically employed in HF and RPA calculations. The 
effective NN interaction which determines the ground-state properties, also deter- 
mines the small amplitude motion around the nuclear ground state. This property 
of the present model ensures that RPA amplitudes do not contain spurious compo- 
nents associated with the center-of-mass translational motion. Models that are not 
fully self-consistent necessitate the inclusion of an additional free parameter in the 
residual interaction, to adjust a proper separation of the spurious state. 

One of the interesting questions is to which extent the UCOM-RPA transition 
spectra are sensitive to the range of the tensor correlator employed in the unitary 
transformation. We have verified that the multipole strength distributions do not 
depend on variations of the central correlator range, around the standard short- 
range correlator 15 . 

In Fig. 1, we display the UCOM-RPA spectra of isoscalar giant monopole 
resonances (ISGMR) for several closed-shell nuclei, using the correlated Argonne 
V18 interaction with different constrains on the range of the tensor correlator, 
I^ S ~ 1,T ~ ^ '=0.07 ', 0.08, and 0.09 fm 3 . For comparison, we also denote excitation en- 
ergies from a selection of experimental 29 > 30 > 31 and theoretical 24 ' 32 ' 33 studies. The 
agreement with experiment and other theoretical results is rather good for the stan- 
dard correlator set with i^' S-1 ' T-0 ' ) =0.09 fm 3 . In heavy nuclei, the ISGMR energies 
are overestimated by « 1 — 3 MeV. By varying the range of the tensor correlator 
around its standard value, the transition strength can be fine-tuned to improve the 
agreement with the experimental data. 

Next we employ the UCOM-RPA to describe the isovector giant dipole reso- 
nances (IVGDR) in 90 Zr, 132 Sn, and 208 Pb (Fig. 2). The correlated Argonne V18 
interaction is used, with different constraints on the ranges of the tensor part of the 
correlator, J^ S_1 ' T_0 ^=0.07, 0.08, and 0.09 fm 3 . The calculated dipole response is 
compared with the experimental data 34 < 35 > 36 > 37 an d with the theoretical excitation 
energies from the relativistic RPA 38 based on DD-ME2 interaction 39 . In all nu- 
clei under consideration, the resulting IVGDR strength distributions display rather 
wide resonance-like structures. The decrease in the range of the tensor correlator, 

/Or rj-i q"i 

i.e. its constraint I# ~ ' ~~ —0.09 fm 3 towards 0.07 fm 3 , results in lower IVGDR 
peak energies by «2-3 MeV. However, UCOM-RPA overestimates the IVGDR cen- 
troid energies by ~3-7 MeV. This difference can serve as a direct measure of the 
missing correlations and three-body contributions in the UCOM-RPA scheme. In- 
clusion of the three-body interaction and long-range correlations beyond the simple 
RPA method, would probably to a large extent resolve the presently obtained dis- 
crepancies with the other studies. 
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Fig. 1. The calculated UCOM-RPA strength distribution of ISGMR for the correlated Argonnc 



(S=l,T=o) 



0.07, 



V18 interaction, using different restrictions on the range of the tensor correlator (7 
0.08, and 0.09 fm 3 ). The experimental data 29,30,31 an( j resm t s f r0 m the nonrelativistic (Drozdz 
et al.) 32 and relativistic RPA (NL3) 24 ' 33 are denoted by arrows. 



In Fig. 3 we show the UCOM-RPA isoscalar quadrupole transition strength dis- 
tributions for 40 Ca, 90 Zr, and 208 Pb (Argonne V18, 4 S=1 ' T=0) =0.07, 0.08, and 0.09 
fm 3 ), in comparison with experimental data 40 . The residual interaction constructed 
from the correlated realistic NN interaction is attractive in the isoscalar channel, 
generating strongly collective peaks corresponding to isoscalar giant quadrupole res- 
onance (ISGQR). In addition, in the case of 90 Zr, and 208 Pb, the UCOM-RPA model 
also results with pronounced low-lying 2 + states. However, RPA based on the cor- 
related realistic NN interaction, without the long-range correlations and three-body 
contributions, is not sufficient for a quantitative description of ISGQR excitation 
energy. For the short range tensor correlator (J^ ~ ' ~ —0.07 fm 3 ) the model still 
overestimates experimental values by w 8 MeV. By decreasing the range of the 
tensor correlator, the quadrupole response is systematically shifted towards lower 
energies. The quadrupole response is rather sensitive to the range of the tensor cor- 
relator. For °Ca and the ranges of tensor correlator determined by 1$ ~ ' ~ =0.07, 
0.08, and 0.09 fm 3 , the ISGQR centroid energies read 25.1, 26.2, and 27.1 MeV, 
respectively. In the cases of heavier nuclei, these differences are smaller, e.g. for 
208 Pb, the centroid energy lowers by 1.2 MeV when going from the correlator with 
7 ( S=1) r=o) =() og fm 3 tQwards 7^ =1 ' T =°) = o.07 fm 3 . 

The agreement achieved between the calculated and experimental properties 
of the ISGMR indicates that the correlated NN interaction corresponds to realis- 
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Fig. 2. The UCOM-RPA strength distributions for the IVGDR in 90 Zr, 132 Sn, and 208 Pb. The 
calculations are based on the correlated Argonne V18 interaction, using different constraints on the 
tensor correlator range (i^ S=1 ' T=0) = 0.07, 0.08, and 0.09 fm 3 ). The experimental data 34,35,36,37 
and the relativistic RPA (DD-ME2) energies 38 > 39 are shown by arrows. 



tic values of the nuclear matter (NM) incompressibility. It has been demonstrated 
in the past that, within relativistic and non-relativistic RPA, the energies of the 
dipole and quadrupole resonances, on one hand, and the value of the effective mass 
corresponding to the effective interaction used, on the other, are correlated 41 ,42 . 
In particular, the relativistic RPA without density-dependent interaction terms, 
based on the ground state with a small effective mass and relatively high com- 
pression modulus, resulted in systematically overestimated energies of giant reso- 
nances 41 . The discrepancies between UCOM-RPA calculations and experimental 
data for multipole giant resonances, as well as the low density of single-nucleon 
UCOM-HF states, suggest that the respective effective mass is too small. Tensor 
correlations with shorter range increase the single-particle level density and result 
in a systematic shift of the giant resonances towards lower energies, improving the 
agreement with experimental data. However, the variations of the ranges of corre- 
lation functions can serve only as a tool for fine tuning of the excitation spectra 
and they can not supplement the effects of the missing long-range correlations and 
three-body contributions. 
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Fig. 3. The ISGQR strength distributions for 40 Ca, 90 Zr, and 208 Pb. The UCOM-RPA model is 

(S=l T=0) 

based on the correlated Argonne V18 interaction with different values of 7^ ' =0.07, 0.08, 
and 0.09 fm 3 to constrain the range of the tensor correlator. The experimental ISGQR excitation 
energies are denoted by arrows 40 . 



4. Conclusions 

In the present study, the fully self-consistent RPA is formulated in the single-nucleon 
Hartree-Fock basis by using correlated realistic NN interactions. The short-range 
central and tensor correlations induced by the NN potential arc explicitly treated 
within the UCOM framework. It is shown that the correlated NN interactions are 
successful in generating collective excitation modes, but for an accurate description 
of experimental data on excitation energies and transition strengths, one has to 
account for the contributions missing in the present treatment. These are (i) long- 
range correlations beyond simple RPA, which can be included within a RPA scheme 
built on the correlated ground state or by including complex configurations within 
Second-RPA, and (ii) induced and genuine three-body interactions, which one could 
try to model by a simple effective three-body force. 
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